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Abstract 

Using stellar population synthesis techniques, we explore the photometric signatures of white dwarf 
progenitor dominated galactic halos, in order to constrain the fraction of halo mass that may be locked-up 
in white dwarf stellar remnants. We first construct a 10 9 Mq stellar halo using the canonical Salpeter initial 
stellar mass distribution, and then allow for an additional component of low- and intermediate-mass stars, 
which ultimately give rise to white dwarf remnants. Microlensing observations towards the Large Magellanic 
Cloud, coupled with several ground-based proper motion surveys, have led to claims that in excess of 20% 
of the dynamical mass of the halo (10 12 M Q ) might be found in white dwarfs. Our results indicate that (1) 
even if only 1% of the dynamical mass of the dark halo today could be attributed to white dwarfs, their main 
sequence progenitors at high redshift (z~3) would have resulted in halos more than 100 times more luminous 
than those expected from conventional initial mass functions alone, and (2) any putative halo white dwarf 
progenitor dominated initial mass function component, regardless of its dynamical importance, would be 
virtually impossible to detect at the present-day, due to its extremely faint surface brightness. 
Keywords: galaxies: halo — galaxies: luminosity function, mass function — galaxies: stellar 
content — stars: white dwarfs 



1 Introduction 

The flat rotation curves of spiral galaxies suggest that their 
galactic halos are mostly composed of dark matter (Rubin, 
Ford & Thonnard 1980). The precise nature of this dark 
matter remains an open question, one with important impli- 
cations for cosmology and galaxy formation. Observations 
of microlensing events from programs such as MACHO and 
EROS suggest that perhaps of order 20% of the Galactic 
halo may be composed of ~0.5 Mq compact objects (Al- 
cock et al. 2000; Afonso et al. 2003). One particularly 
appealing source of 0.5 Mq objects that was immediately 
vetted in the literature was that of a population of faint, 
old, white dwarfs (WDs). 1 

Alcock et al. (2000) analysed 5.7 years of photomet- 
ric data on stars in the Large Magellanic Cloud, in search 
of gravitational microlensing events. They concluded that 
the fraction of massive compact halo objects (MACHOs) 
is about 20% (with an average mass of 0.5 Mq per MA- 
CHO), and that the total mass in MACHOs out to 50 kpc 
is 10 11 Mq , assuming the lenses are located in the Galactic 
halo. 2 The recent EROS data toward the Small Magellanic 
Cloud allow for a maximum of ~25% of a spherical, isother- 
mal, and isotropic Galactic halo of mass 4xlO n Mq out to 
50 kpc which could be composed of objects with mass be- 
tween 2xl0~ 7 M and 1 Mq (Afonso et al. 2003). The ab- 
sence of events with crossing times shorter than 10 days and 
the lack of sufficient numbers of low-mass main sequence 
stars in the Hubble Deep Field essentially rules out planet- 

1 "Appealing" in the sense that they are relatively "mundane" 
astronomical objects, obviating the need to draw upon more ex- 
otic solutions. 

2 Not necessarily a unanimously accepted hypothesis - e.g. 
Sahu & Sahu (1998); Di Stefano (2000). 



like objects and brown dwarfs as the microlens candidates 
(Gould, Bahcall & Flynn 1998; Lasserre et al. 2000). The 
current favoured mass range from the microlensing exper- 
iments also rules out the remnants of more massive stars, 
including neutron stars and black holes, as potential (and 
substantial) baryonic dark matter candidates. 

A rich literature has emerged over the past five years 
exploring the pros and cons of the hypothesis that white 
dwarfs may comprise a significant fraction of the dynamical 
mass of the galaxies. Ryu, Olive & Silk (1990) constructed 
a simple galactic halo chemical evolution model and ruled 
out neutron stars, but not white dwarfs, as dark matter 
candidates, based on metallicity and luminosity considera- 
tions. The low incidence of Type la supernovae observed in 
the outer regions of galaxies lead Smecker & Wyse (1991) to 
impose tight upper limits on the fraction of white dwarfs in 
galactic halos. Chariot & Silk (1995), using population syn- 
thesis models coupled with number counts of faint galaxies 
in deep galaxy surveys, showed that only a small fraction 
(<10%) of present-day halos could be in the form of white 
dwarfs. Madau & Pozzetti (2000) showed that in order to 
avoid the overproduction of extragalactic background light, 
the halo white dwarf mass fraction should be <5%. Gib- 
son & Mould (1997) and Brook, Kawata & Gibson (2003) 
placed more severe constraints upon the halo white dwarf 
mass fraction (<l-2%) by using the observed carbon, ni- 
trogen, and oxygen abundance patterns of halo stars. 

Following Larson's (1986) suggestion that a remnant- 
dominated form of the initial mass function (IMF) could 
account for the unseen mass in the solar neighborhood, both 
Chabrier, Segretain & Mera (1996) and Adams & Laugh- 
lin (1996) devised physically- motivated models allowing for 
a white dwarf progenitor dominated IMF, consistent with 
the aforementioned MACHO and EROS microlensing re- 
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suits. Such an IMF differs from that of the classical Salpeter 
(1955) functional form through the absence of both low 
mass (<1 M ) and high mass (>6 M ) stars (Gibson & 
Mould 1997; Fig 1). 

In this paper, we address explicitly the photometric 
properties of the luminous early phases of putative white 
dwarf progenitor dominated galactic halos, in addition to 
their temporal evolution. In Section 2, we describe our stel- 
lar population models and the different functional forms for 
the IMFs considered here. The results of our calculations 
are then presented in Section 3, with the accompanying 
discussion and conclusions provided in Section 4. 



2 Models 

Using our evolutionary stellar population synthesis code 
(Lee, Lee & Gibson 2002), we have calculated the evolu- 
tion of the photometric properties of galactic halos as a 
function of the dynamical mass fraction tied up in putative 
populations of white dwarfs. We allow the total present- 
day mass fraction of white dwarfs to range from - 100% 

3 of the halo's dynamical mass, which is taken for this ex- 
ercise to be 10 12 M (after that of the Milky Way - Fich 
& Tremaine 1991). All models described herein contain 
a standard 10 9 M stellar component 4 , which is itself de- 
scribed by the Salpeter IMF (by number): 

$(m)dm = dn/dm — AmT x dm, (1) 

with £=2.35. 

The fraction of halo dark matter contained in white 
dwarfs is then determined by the amount of matter con- 
tained in a supplemental IMF described by a truncated 
power-law of the form 

<E>(m) dm = dn/dm = Ae'^" 1 ^ x m~ Q dm (2) 

for which we use m = 2.7, (3 = 2.2, and a = 5.75 (Chabrier 
et al. 1996). The peak of this skewed-Gaussian functional 
form for the IMF (hereafter, wdlMF) occurs at m f«2 M Q , 
favouring the production of white dwarf progenitors. The 
wdlMF yields a present day Galactic halo mass-to-light ra- 
tio S> 100 after a Hubble time as most of its initial stellar 
distribution has since become very faint remnants. 

The stellar population synthesis models presented here 
are based upon the Y 2 Isochrones 5 (Kim et al. 2002) with 
[a/Fe]=+0.3, coupled to the post-red giant branch stellar 
evolutionary tracks of Yi, Demarque & Kim (1997). We 
have calibrated the horizontal-branch morphology with the 
Milky Way globular clusters, as in Lee et al. (2000, 2002). 
The stellar library of Lejeune, Cuisinier & Buser (1998) 
was taken for the conversion from theoretical to observable 
quantities. In our calculations, the stellar remnant mass 
was assumed to be 0.5 M for initial stellar masses below 
8 M , and 1.0 M Q for initial stellar masses in excess of 
8 M . 

3 Strictly speaking, we allow up to 99.9% present-day mass 
fraction of white dwarfs as a 10 9 M stellar halo with the canon- 
ical Salpeter IMF provides the base for our model. 

4 A stellar mass of 10 9 M for the Milky Way halo is rea- 
sonable once one adopts, for example, a surface density profile 
proportional to r~ 2 normalized by the local stellar halo density 
- e.g. Preston ct al. (1991). 

5 http : //csaweb.yonsei . ac .kr/~kim/yyiso .html 



New cooling models for white dwarfs have become avail- 
able recently which possess somewhat different behaviour in 
the colour-magnitude plane at advanced ages (e.g. Hansen 
1998; Richer et al. 2000), when compared with the conven- 
tional models adopted in most population sysnthesis codes 
(including ours). While we have experimented with their in- 
clusion, due to the extremely low luminosity at which these 
new models diverge from the classical ones (at Mv~17.5, 
after ~8 Gyr of cooling), our results are not seriously im- 
pacted by the choice of specific white dwarf cooling tracks. 

The following section presents results for our default 
model - a Milky Way-like halo with metallicity [Fe/H] ~ 
— 1.8 (Ryan & Norris 1991) - we also touch briefly upon 
the implications for more metal-rich halos such as that for 
M31. 



3 Results 

Figure 1 shows the variation in absolute V-band magni- 
tude as a function of halo dynamical mass fraction tied 
up in white dwarfs, for ages ranging from 1 to 14 Gyr. 
The asymptotic values of My approached for mass fractions 
<10~ 4 simply reflect the model's assumed base 10 9 M 
Salpeter IMF stellar component; 6 a mass fraction of unity 
corresponds to a 10 12 M halo comprised entirely (essen- 
tially) of white dwarfs. After ~13 Gyr, a 100% white 
dwarf halo is, for all intents and purposes, indistinguishable 
photometrically from a halo containing no remnants aside 
from the underlying population attributed to the (known) 
Salpeter-like component. For younger stellar systems though, 
WD progenitors preferentially populate the main sequence, 
giant branch, and horizontal branch phases and rapidly 
dominate the luminosity of the halo. 

Although unlikely, if the total dynamic mass were com- 
posed of white dwarf remnants, this halo would have been 
~10 magnitudes brighter ~1 Gyr after its "formation" (i.e. 
z~3, for a redshift of formation of five). In comparison, 
the same halo (by dynamical mass), but now with only 
the known 10 M Salpeter IMF stellar component, would 
only be ~2 magnitudes brighter at the same lookback time. 
Even if only 1% of the mass of the dark halo was in the form 
of this wdlMF component, it would have been a factor of 
~100 times more luminous than the canonical halo at the 
same redshift z~3. 

Figure 2 presents the dependence of colour (V— I) on 
the fraction of the dynamic halo mass in the form WD 
remnants for halos ranging in age from 1 to 14 Gyr. A 
white dwarf dominated halo is indistinguishable in colour 

6 We re-emphasie here that the asymptotic value of My ap- 
proached at low white dwarf mass fractions in Figure 1 corre- 
sponds to the values expected for a halo with a base stellar 
component of mass 10 9 M , constructed with a conventional 
Salpeter IMF - representative of what one might expect for a 
canonical model of the Milky Way. This absolute magnitude 
would of course change should a different "base" Salpeter IMF 
halo component be employed - for example M.y for a halo of the 
same metallicity would be 2.5 mag brighter should the base stel- 
lar halo component have a mass of 10 10 M . Furthermore, the 
halo would be about 1 mag fainter if its metallicity was solar, as 
the main-sequence turnoff is fainter for metal-rich systems. Re- 
gardless, our analysis remains valid, as it is the differential be- 
tween the models which is important, and not the "zero-point" 
of the curves in Figure 1. 
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Figure 1: The absolute V-band magnitude as a func- 
tion of halo white dwarf mass fraction for eight ages 
ranging from 1 to 14 Gyr. The left-hand limit to the 
plot corresponds to the minimum contribution due to 
the known 10 9 M Q stellar halo intrinsic to spirals such 
as the Milky Way (our default model). Passive evo- 
lution of this basic system would be ~2 magnitudes 
brighter at redshift (an age of ~1 Gyr, for a red- 
shift of formation of five). Conversely, a 10 12 M Q halo 
comprised of nothing but white dwarfs at the present- 
day would have been ~10 magnitudes brighter at z~3. 
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Figure 2: As in Figure 1, but now showing the variation 
in integrated colour (V— I) as a function of halo white 
dwarf mass fraction. 



from a Salpeter-only halo for lookback times <8 Gyr (note 
the essentially flat behaviour of the curves as a function of 
wdlMF mass fraction). 7 

The only marginal difference between the two scenar- 
ios occurs at ages corresponding to redshifts in the range 
l<z<2 - for instance, if 1% of the dark halo mass was in the 
form of WDs, its stellar component would have been ~0.08 
mag redder than a Salpeter-only halo (left-hand limit of 
Figure 2) at early epochs because of the former's more sig- 
nificant populations of red giants and red horizontal-branch 
stars. 

The predicted colours in Figure 2 are also, of course, 
different for different metallicities. In general, the more 
metal-rich a halo is, the redder its stellar population. For 
a halo with a metallicity like that of M31 ([Fe/H] ~ -0.5 
- Durrell, Harris & Pritchet 2001), we note in passing that 
if 1% of the halo's dynamical mass were attributed to a 
wdlMF component, the V— I colour at 1 Gyr would be 0.55, 
at 3 Gyr it would be 1.14, and at 14 Gyr it would be 1.26. 



4 Discussion 

In Figure 3 (after Chabrier 2001), we speculate on the sur- 
face brightness radial profile of a 10 12 M0 halo with a 
10% (by mass) wdlMF component, but without the "base" 
10 9 M Salpeter IMF stellar component. 8 We have adopted 
a dark matter halo density profile of p(r) cx 1/r 2 . For the 
surface mass density of the halo, we assumed 23 Mq/pc 2 
at 8 kpc (Kuijken & Gilmore 1991). It can be seen from 
Figure 3 that if 10% of the dark halo mass were composed 
of a wdlMF stellar component, its predicted surface bright- 
ness would be essentially undetectable. Put into context, 
the central surface brightness of the low surface brightness 
galaxy Malin 1 is ^ (V)~25.5 (Bothun et al. 1987) - i.e., 
this 10% wdlMF component lies 7-10 magnitudes fainter 
than the surface brightness of Malin 1! 



Table 1: Surface Brightness at 60 kpc (Z=0.0004) 



Age (wdlMF) 


B 


V 


R 


I 


12 Gyr (1%) 


38.70 


37.44 


36.59 


35.78 


12 Gyr (10%) 


36.20 


34.94 


34.09 


33.28 


12 Gyr (100%) 


33.70 


32.44 


31.59 


30.78 


14 Gyr (1%) 


40.45 


39.12 


38.22 


37.38 


14 Gyr (10%) 


37.95 


36.62 


35.72 


34.88 


14 Gyr (100%) 


35.45 


34.12 


33.22 


32.38 



Table 1 lists the surface brightness of such a galac- 
tic halo at a galactocentric radius of 60 kpc for several 
passbands at a given age and wdlMF mass fraction. We 
do not wish to belabour the result here, and simply note 
that the surface brightnesses in the outer regions of these 
putative WD-dominated halos would be essentially unde- 
tectable. More metal-rich halos would be somewhat brighter 

7 The abrupt shift to significantly bluer colours for the 1 Gyr 
isochrone is due to the treatment of convective core overshooting 
in the stellar models (Yi 2003; Fig 6). 

8 By eliminating the underlying trace (by mass, not by lu- 
minosity) Salpeter component, we serve to emphasise the light 
contribution from the wdlMF component alone. 
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distance (kpc) 



Figure 3: Surface brightness radial profile of a galactic 
halo with 10% of its dynamical mass comprised of a 
wdlMF stellar component, but now without the base 
10 9 Mq Salpeter IMF stellar component - to empha- 
sise the essentially undetectable wdlMF component's 
surface brightness - after Chabrier (2001). 

because the main-sequence lifetime of stars increases with 
stellar metallicity, but still well below any realistically ob- 
servable limit. 



5 Summary 

We have explored the photometric ramifications of ascrib- 
ing a significant fraction of a galactic halos dynamical mass 
to a population of old, faint, white dwarfs. Using our pop- 
ulation synthesis package, we show that while such an hy- 
pothesis is essentially impossible to rule out using surface 
brightness arguments for nearby (old) halos, passive evo- 
lution of the progenitors of any large, putative, halo white 
dwarf population would make the present-day halos 100 - 
500 times brighter at redshifts 1< z <3 than halos evolving 
passively under standard stellar evolution arguments (i.e. 
a Salpeter- like initial mass function). 

Making one-to-one comparisons between single objects 
such as the Milky Way and potentially quite different ob- 
jects at high-redshift can be dangerous; what we wish to 
end with is simply a note drawing attention to, for ex- 
ample, the existence of passively-evolving L* galaxies at 
high-redshift which individually are consistent with passive 
evolution of normal stellar populations and normal IMFs 
(e.g. Waddington et al. 2002) - i.e. the stellar popula- 
tions therein are < 2 mags brighter then than today (un- 
der reasonable assumptions of mass, etc), opposed to being 
> 4-5 mags brighter than expected if these galaxies were 
dominated dynamically by white dwarf precursors. A safer 
comparison is to model the ensemble of the population (e.g. 



extragalactic background light and deep number counts), 
both of which are also inconsistent with the white dwarf 
scenario (Chariot & Silk 1995; Madau & Pozzetti 2000). 
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